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A  COMPUTER  PROGRAM  AND  APPROXIMATE  SOLUTION  FORMULATION 
FOR  GUN  MOTIONS  ANALYSIS 


ABSTRACT .  Jhe  purpose  of  this  paper  is  to  describe  some  of  the  features 
associated  with  a  finite  element  computer  program  for  approximate  solutions 
of  a  gun  dynamics  problem.  The  lateral  motion  of  a  gun  tube  is  modeled  by 
an  Eul er- Bernoul li  beam.  The  difficulties  of  the  problem  are  due  to  various 
complicated  loadings  and  support  conditions  which  can  be  nonconservative, 
highly  discontinuous  and  time  dependent.  The  solution  formulation  for 
this  generally  non-self-adjoint  problem  has  been  presented  in  an  earlier 
paper.  In  terms  of  finite  element  discretization,  the  two-dimensional 
shape  function  of  spatial  and  time  coordinates  is  chosen  as  a  product  of 
two  one-dimensional  shape  functions;  each  for  its  respective  coordinate 
abd  both  being  Uermitian  polynomials.  The  generalized  coordinates  are 
then  the  displacement,  slope,  velocity  and  time  derivatives  of  the  slope  at 
each  node  point.  The  correspondence  between  local  and  global  generalized 
coordinates  is  described.  The  "stiffness  matrices"  of  spatial  and  time- 
effect,  contributed  by  the  recoil  force,  pressure  and  curvature  induced 
force  and  the  moving  mass  of  a  projectile  are  derived.  It  is  interesting 
to  observe  that  the  strong  discontinuities  associated  with  these  forces 
disappear  as  a  result  of  the  smoothing  effect  of  integration  in  spatial 
as  well  as  in  time  coordinates.  The  present  approach  to  deal  with  the 
moving  support  problem  efficiently  is  also  pointed  out  in  this  paper. 
Numerical  results  of  a  demonstrative  problem  are  presented. 

1 »  I RODUCT I ON .  It  is  our  ultimate  goal  to  design  a  gun  system 
from  which  a  projectile  can  hit  a  target  with  precision.  The  technology 
involved  in  this  task  covers  that  of  the  interior  and  exterior  ballistics. 

As  it  is  obvious  from  the  sequence  of  events,  the  position,  direction  and 
the  velocity  of  a  projectile  as  it  leaves  the  muzzle,  consititue  the 
initial  conditions  for  the  problem  of  exterior  ballistics.  Clearly  then, 
the  capability  of  providing  these  data  under  a  variety  of  loading,  support¬ 
ing  and  detonating  conditions  is  one  of  the  essential  links  for  the 
complete  chain  of  designing  a  precision  gun  system. 

This  paper  reports  some  of  our  results  in  the  initial  effort  on  the 
gun  motions  analysis.  Specifically,  these  are  some  of  the  detailed  features 
related  to  a  computer  program  for  such  an  analysis.  The  finite  element 
method  in  a  very  general  sense  has  been  chosen  here  due  to  its  proven 
flexibility  and  adaptability  for  a  very  wide  range  of  situations  in  terms 
of  geometry,  loading  and  supporting  conditions.  Even  at  this  early  stage. 
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we  have  observed  some  unusual  features  of  the  problem  which  render  the  use 
of  "canned"  finite  element  computer  program  meaningless.  First,  there  is 
no  satisfactory  finite  element  computer  program  for  initial  value  programs 
at  the  present  time  other  than  the  methodology  developed  by  this  writer 
using  an  unconstrained  variational  formulation  [1] .  The  forces  considered 
in  the  present  study  include  the  moving  mass  effect  of  the  projectile, 
the  curvature  and  pressure  induced  load  and  the  recoil  force  which  is 
nonconservative  in  nature  if  the  support  at  the  breech  end  is  less  than 
rigid.  Various  elastic  supports  are  considered.  In  the  due  course, 
effects  of  frictional  forces,  large  deformation  and  time-dependent 
support,  and  other  effects  will  be  included.  Therefore,  it  is  only  sen¬ 
sible  to  develop  right  from  the  beginning  a  finite  element  computer  pro¬ 
gram  codes  uniquely  suitable  for  the  problem  on  hand,  which  is  capable 
of  handling  nonconservative  forces,  moving  mass  effects,  various  support¬ 
ing  conditions  including  moving  supports.  In  addition,  it  should  also  be 
easily  amenable  for  modifications  and  inclusions  of  other  effects  such  as 
frictional  forces,  large  deformation,  and  so  forth. 

The  unconstrained  variational  formulation  for  approximate  solutions 
of  boundary  value  problems,  in  conjunction  with  the  finite  element  descrit- 
ization  has  proved  to  be  highly  efficient,  easy  to  use  and  utterly  flexible 
for  problems  involving  nonconservative  forces,  various  support  conditions 
and  damping  effects  [2,3].  Its  unique  advantage  in  dealing  with  initial 
value  problems  has  also  been  demonstrated  [1] .  Therefore  it  is  natural 
that  this  finite  element  unconstrained  variational  formulation  be  adapted 
for  the  complicated  transient  problem  of  gun  motions  analysis. 

The  basis  of  the  present  formulation  for  more  general  cases  has  been 
given  in  a  previous  paper  [4] .  The  special  problem  of  a  uniform  gun  tube 
is  treated  here  for  demonstrated  purposes.  The  differential,  initial  and 
boundary  conditions  are  given  in  Section  2.  This  differential  equation 
differs  from  those  used  by  many  other  writers  by  one  term.  This  extra 
term  included  in  previous  papers  turned  out  to  be  incorrect  and  the 
reason  is  given  in  Section  3.  An  unconstrained  variational  statement 
which  is  equivalent  to  the  given  governing  equation  is  stated  in  Section  4. 
The  details  of  some  of  the  special  features  on  finite  element  descritization 
are  described  in  Section  5.  Finally  some  preliminary  data  which  shows  the 
recovery  of  the  initial  data  are  presented  in  Section  6. 

2.  GOVERNING  EQUATIONS.  The  motion  of  a  gun  tube  modeled  by  the 
lateral  deflection  of  an  Euler-Bernoulli  beam  is  shown  in  Figure  1.  The 
differential  equation  in  nondimensional  form  is 

y""  +  (-  P  +  g  sina) [ (l-x)y ' ] '  +  y 


=  -  Py"  H(x-x) 

-  m[g2t2y"  +  2gty'  +  y]6(x-x) 
-  (gm  cosa)6(x-x)  -  g  cosa 


(2-1) 


with 


P 


=  irR2p 


2 


i 


FIGURE  1.  A  Schematic  Drawing  of  the  Problem  Configuration. 
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where 


y  =  y(x>t),  the  tube  deflection 

x,  spatial  axis  along  the  tube's  length,  0<x<l 

t,  time  axis,  0<t<T,  T  is  the  time  limit  of  Tnterest 
a,  elevation  angle 
m,  projectile  mass 

3,  acceleration  of  the  projectile,  assumed  to  be  constant 
p,  bore  pressure,  assumed  to  be  constant 
g,  gravitational  acceleration 

-  _  1  2 

x  -  2  Bt  ,  projectile  position 

H(x),  the  Heavenside  step  function 
6(x),  the  Dirac  delta  function 


"n  j  L*  ’  a.Pr:*-me  O  denotes  a  differentation  with  respect  to  x  and  dot  (*), 
a  ^fferentiatron  with  respect  to  t.  The  derivation  of  this  equation  and 
ie  end  conditions  which  follow  have  been  given  previously  [4,5]  and  will 
not  be  repeated  here.  However,  the  absence  of  a  term  in  the  present  equa¬ 
tion  will  be  explained  in  the  next  Section. 


The  initial  condition, 
time  are 


or  more  appropriately,  the  end  conditions  in 


y(x,0)  =  0 


(2-2) 


y(x,T) [1  +  mfi (I  3T2-x) ]  +  k?[yCx,0)  -  Y(x)]  =  0  . 

where  the  constant  k7  is  introduced  in  conjunction  with  the  unconstrained 
variational  formulation  (Section  4)  so  that  if  one  takes  k7  to  be  infinite 

dlsPlacemen^  y(x,0)  is  forced  to  be  identical  to  the  prescribed 
s  dpe  (x) .  based  on  similar  reasonings  the  boundary  conditions  have  been 
shown  to  be  the  following. 


and 


y"(0,t)  -  k2y*(o,t)  =  0 
y”(l,t)  +  k4y'(l,t)  =  0 


y'"  (0,t)  +  kxy (0,t) 

+  (-  P  +  g  cosa)y'(0,t)  +  Py' (0,t)H(I  Bt2) 
+  mg2y' (0,t)5(j  Bt2)  =  0 
y"'  (l>t)  -  k3(l,t)  +  Py'(l,t)H(i  6t2-l) 

+  mB2y' (l,t )6(i  Bt2-1)  =  0 


(2-3) 
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I 

where  ki}  i  =  1,2, 3,4,  are  the  appropriate  elastic  spring  constants  at 
the  supports. 

3,  AN  ERROR  IN  PREVIOUS  LITERATURE*.  We  have  already  mentioned 
that  there  should  be  a  correction  made  to  the  force  terms  due  to  the 
projectile's  moving  mass.  This  correction  will  be  explained  in  the 
present  section. 


The  geometry  of  a  concentrated  mass  moving  on  a  "beam"  structure  is 
illustrated  in  Figure  2.  In  this  figure,  y(x,t)  is  the  beam  deflection. 

The  position  vector  r  of  the  moving  mass  m  can  be  written  as 

r  =  £(t)i  +  y(£(t)  ,t)  j  (3-1) 

where  i,j  are  the  unit  vectors  in  s  and  y  directions  respectively;  £(t) 
denotes  the  position  of  m  in  x  direction.  It  is  easy  to  derive  the  velocity 

and  acceleration  of  r. 

H.  *> 

i  = 'd=r-  £(t)i  +  [C(t)y'(5(t),t)  +  SL  ]j 

l  *  dt  =  C(t)i  +  KCOy'CC.t)  +  S2(t)y"(S,t)  I  (3_2) 

+2|(t)  ]j 

3£3t  3t2  - 

In  deriving  (3-2),  it  should  be  observed  that 


4-y(C(t),t) 


=  &  +  9y 

3£  3t  3t 


(3-3) 


’‘This  error  was  first  observed  by  T.  E.  Simkins  of  Benet  Weapons  Laboratory 
at  the  Second  U.S.  Army  Symposium  on  Gun  Dynamics,  19-22  September  1978, 
Rensselaerville,  New  York. 
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Since  a  is  the  acceleration  of  the  projectile  mass  m,  the  force  acting  on 
it  must  be 


and 


F  =  ma 


Fx  =  max  = 


=  m{£y' 


+  £2y"  +2^y< 


(3-4) 


It  appears  that  most  of  the  previous  investigators  have  accomplished  this 
much.  Now  comes  the  important  step,  since  we  are  interested  in  the  force 
imparted  on  the  beam  by  the  moving  mass  m.  This  force,  from  Figure  2,  is 
clearly  not  due  to  Fy  alone,  rather  is  the  sum  of  the  components  of  Fx  and 
Fy  in  the  normal  direction  of  the  beam.  Denote  this  sum  by  Fm,  one  has 


Fm  =  Fy  cos0  -  Fx  sin0 

But,  for  small  deflection,  0  is  small  so  that 

cos0  5  l 

sin0  =  tan0  =  0  =  =  y ' 

35 


(3-5) 


(3-6) 


Substitute  (3-6)  into  (3-5),  one  obtains 

Fm  =  F  -  m£y* 
or 

Fm  =m^2y"  +25 

at  at2 

Eq.  (3-7)  is  the  correct  expression  for  the  force  imparted  on  a  carrying 
structure  due  to  a  moving  mass.  This  structure  can  be  a  rail,  a  bridge 
or,  as  in  the  present  case, a  gun  tube.  It  is  somewhat  surprising  to  find 
that  the  expression  for  Fy  of  Eq.  (3-4)  rather  than  Fm  of  Eq.  (3-7)  is 
used  in  the  literature  known  to  this  writer  (see,  for  example,  [6]). 


r 

(3-7) 


6 


«t) 
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FIGURE  2,  Forces  Acting  Between  Gun  Tube  and  a  Moving  Projectile. 


4.  UNCONSTRAINED  VARIATIONAL  STATEMENT.  Through  integrations-by- 
parts,  it  is  an  easy  matter  to  show  that  the  following  variational  state¬ 
ment  is  equivalent  to  the  differential  equation  and  ends  conditions  stated 
in  Section  2. 


with 


and 
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61  =  (61)  =  l  C6I±)  -  l  (6J • )  =  0 

i=l  7  j=l  J 

1  1 

J,  =  -  Tg  cos af  J  y*dxdt 
0  0 

1  1  _  _ 

^2  ~  ~  cosa /  /  y*6(x-x)dxdt 

0  0 

1 

J3  =  kyf  Y(x)y*(x,l)dx 


1  1 

II  =  /  /  y"y*"dxdt 

0  0 

1  1 

I2  =  (p  -  §  sina )/  /  (l-x)y'y*'dxdt 

0  0 

1  1.. 

I,  =  -  /  /  yy*dxdt 

0  0 

-11 

I.  =  -  Pj  j  y'y*'  H(x-x)dxdt 
0  0 

_  1  1 

I_  =  -  pj  /  y’y* '6(x-x)dxdt 
b  0  0 


(4-1) 


(4-2) 


^(4-3, first  part) 
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Q2  1  1  -  _  \ 

I  =  -  !B| -J  /  tVy*'<5(x-x)dxdt 
0  0 

«2  1  1 

I  =  -  |  {  t  y'y* ' 6 1 (x-x)dxdt 

7  TOO 
-  R  1  1  . 

Ig  =  t  /  /  ty'y*6(x-x)dxdt 
1  0  0 

a  1  1 

Ig  =  ?■  /  J  y'y*<5(x-x)dxdt  l 

y  1  0  0 

1  1..  -  - 

I  =  -  EL  J  J  yy*6(x-x)dxdt  (4-3, second  part) 

10  T  0  0 

li.:. 

I„--y/  /  yy*S(x-x)dxdt 

1  1  0  0 

1 

i12  =  +  k2y,(0,t)y*,(0,t) 

+  k3y(l,t)y*(l,t)  +  k4y'(l,t)y*'(l,t) 

+  k5y(xs,t)y*(xs,t)  +  k6y'(xs,t)y*'(xs,t)}dt 

1 

Il3  =  k7/  y(x,0)y*(x,l)dx  / 

0 

This  variational  statement  will  serve  as  the  basis  of  our  finite  ele¬ 
ment  solutions. 

5.  FINITE  ELEMENT  DISCRETIZATION. 

I 

5.1  From  Local  to  Global  Coordinates.  The  purpose  of  the  discretiza¬ 
tion  is  to  enable  one  to  write  the  variational  statement  of  Eq.  (4-1),  which 
is  a  functional  of  continuous  functions  y  and  y*,  etc.,  in  the  form  of  a 
matrix  equation  < 

6Y*T  K  Y  =  6Y*T  F 


(5-1) 


(l+JOCl+'l)1 - 1 - - - 1 - - 1 - 1 - - - 1 - -  i+Ci+m)! 
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FIGURE  3,  Globe  Node  Scheme  for  the  (ra,n)th  Element. 


where  Y,  Y*  are  the  "global"  generalized  coordinates  vectors.  K  is  the 
global~"stiffness"  matrix,  and  F  the  "force"  vector.  These  terminology 
are  patented  after  the  static  structural  analysis,  but  they  do  not  nec¬ 
essarily  have  the  physical  meanings  of  those  adjectives  attached  to  them. 
Since  the  variational  statement  associated  with  (5-1)  is  unconstrained,  the 
equation  leads  directly  to 


(5-2) 


K  Y  =  F 


which  can  be  solved  for  Y  if  K  and  F  are  properly  defined.  The  process  by 
which  K  and  F  are  assembled  and  the~relation  between  Y  and  the  desired  solu 
tion  y”(x,t)  will  be  described  in  detail  here  in  this  section. 

The  first  step  is  to  write  down  the  expressions  in  the  variational 
statement  in  terms  of  the  element  variables.  A  grid  scheme  of  elements  is 
shown  in  Figure  3.  In  this  figure,  the  nondimensional  length  of  the  gun 
tube  is  divided  into  K  equal  segments  and  the  time  range  of  interest  into 
L  equal  segments.  The  result  is  then  a  set  of  K  x  L  rectangular  elements. 
In  the  equations  that  follow  n,  the  sub-  or  super-scripts  m,n  denote  the 
association  with  the  m^1,  n^h  segments  or  the  (m,n)th  element.  Define  the 
relation  between  the  local  coordinates  (£»n)  of  the  (mjn)^  element  and  the 
global  coordinates  (x,t)  by 


£  =  £(m)  =  Kx  -  m  +  1 
q  =  qCn)  =  Lt  -  n  +  1 


(5-3) 


Or 


x  =  I  (S+m-1) 
K 


(5-4) 


t  =  —  (q+n-1) 
L 


One  can  write,  from  Eqs.  (4-2)  and  (4-3): 


(S-5a) 


11 


1  1 


'N 


{I3  -  -  Sif  lo  >>, „)<*(„, n)^dn 


{i4  =  - St<  s:.,)H(».n)«-«d5dn 


1  1 


1  1 


"s-nsp/  »)W) 


«6  -  -  /1[Cm-i)!*2(u-i)n*nJ])r(B  n)5y(^  n)6(ln  j  (j-5)d5dn 

L  0  0 


32^3^2  1  1 


6i7  =  -  H  mg  1  —  /  /  [(m-D^Cm-Dn+nMyJn^jfiyf^njfifn^^K-OdWn 


0  0 


«8  *  II 


2mgTK 


1  1 


t  I0  (,IC"-1)*'ll>’(n..n)S>'Cn,,n)'5Cm,n)«-OdS<in 


«9  *  II  =?*  /  /  >"(m,n)5ykn){(m,n)<5-0«dn 


1  1 


0  0 


1  1 


SI10  *  -  Hr  4  f/(m.n)S*(m,n)S(m,n)«-»dSdn 


61,-  -  -  XI  —  /  I  Y,  -Ar  A  ,(C-OdCdn 
11  Too  (m>n)  (m»n)  (m»n) 


1  1 


'll 


{I12  •  IfilVo,,)1';.,,)  *  V^U.nlSU 


VtK.nJ^CK.n)  *  k4KV,n)5y[;K,n)ldn 

“H  *  Ji  T  ^,Ds4,l)dc 


C5- 


/ 


5b) 
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and 


K  L  _  1  1  . 

■1  -  -  t  L  ^  I  I  *  n,  »■>*** 
111=1  11=1  KL  0  0  1  •  1 


« i  =  -  Z  I 


SJ3  =  J1TiY00«,^.L)d« 


(5-6) 


Now,  the  shape  functions  are  introduced.  Let 

y(m,n)(?’n)  = 

y(m,n)^,n-)  =  aT(?,n)Y*fm»n) 


(5-7) 


where  a(E,ri)  is  the  shape  function  vector  to  be  chosen  and  Y(m>n)  is  the 
generalized  coordinates  vector  associated  with  the  (m.n)^  element.  The 
choice  of  a(£,ri)  and  the  meaning  of  Y(m>n)  will  now  be  given.  We  shall 
further  write 


re.o)  - 

>V„)  «■">  ■ 

Here  in  this  paper,  we  select  ak(£,r|)  as 

ak(5,n)  =  bi(^)bj (n) , 

where  b^(^)  or  bi(p)  are  defined  by 


l  ak(C,n)Yk(m’n) 

k=l 

k=l  K  K 

i,j=l, 2,3,4 
k=l , 2, ...  16 


bx (O  =  1  -  +  2£3 

b2co  «  k  -  2?2  +  e 

b3(C)  =  3£2  -  2£3 

b4CO  «  -  K2  +  Z3 


(5-8) 


(5-9) 


(5-10) 
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Node  3 


Node  1 


FIGURE  4.  Element  Node  Scheme. 
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TABLE  I.  CORRESPONDENCE  BETWEEN  LOCAL  AND  GLOBAL  NODE  NUMBERS 


Node  Number 

Local 

Global 

1 

Cn-1) (K+l)+m 

2 

(n-1) (K+l) +m+l 

3 

n(K+l)+m 

4 

n(K+l)+m+l 

Now,  the  local  node  scheme  is  defined  in  Figure  4  and  the  correspondence 
between  a  local  node  of  the  (m,n)^  element  and  a  global  node  is  given  in 
Table  I. 

With  Eqs.  (5-9)  and  (5-10)  in  conjunction  of  the  equation  of  discretiza¬ 
tion  (5-8),  the  correspondence  between  the  index  k  and  the  pair  (i,j)  of  equa¬ 
tion  (5-9)  and  that  between  (m,n)  and  y(m,n) (S,0)  can  be  easily  established. 
These  relations  are  given  in  Table  II. 
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In  terms  of  a(£,r)),  Y  and  Y*,  the  expressions  of  61.,  6J.,  etc.  can  be 

written  as  11 


61 


£  k  t,,3  t  11 

1  =  I  I  ~T~  6Yfm  nl  /  /  a  FFa  FF  d^dn  Y 

m=l  n=l  L  0  0 


~(m,n) 


K  L  _  t  11 

617  =  I  I  (P-gsina)  f  6Y*m  ^  /  /  [K-m+l-£ja  d£dn  Y 


m=l  n=l 


L  ~(m,n)  JQ 


(m,n) 


K  L  11 

51 3  =  I  I  ("  -■)  /  /  a  aT  d£dp  Y  (  . 

n=l  n=l  TK  o  0  ~,TV”,T1  ~(n>,n) 


m= 


SI4  -  n  (-  HC.n)«-«<*5dn  v(n>n) 

«I5  -  II  (  (  i  i*  \n,n)d-aw n  Ytm>n) 


5I6  *  II  <-  2&T3£i)«Y*T  /  /*  [tm-l)!.2(m-l)n.ni] 


(m,n) 


0  0 


«I7  «  II  c- 


a  aTr  6.  (C-C)dCdn  Y,  , 
(m,n)  w  ^  ~(m,n) 


1  1 


g£.  T  K  -)6Y*^  >  f  {  [(m-l)2+2(m-l)n+n2]  • 

l3  ~im,nj  Q  Q 

•  a  aTr  6*  (l-C)dCdrj  Y  , 

~  ~*S  (m,n)  ^  *  ~(m,n) 


1  1 


61  =  H  (2mg  ^)6Y**  /  /  [(a-l)+n]a  a  £  6  (£-S)d£dn  Y 

(m,nj  o  0  ^  (m,n)  ‘“(m.n) 


(5-11, part  A) 
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61  =  yy  n*§2JL  6Y*‘  .  /  /  a  a1.  6  ,  (5-5)d5dn  Y  - 

9  LL  L  ~(rn,n)  JQ  Jq - »?  (m,n)  ~(m,n) 


1  1 


61  =  Yy  (-  !0k)5Y*‘  /  /  a  a*  6  (£-£)d5dri  Y 

10  t  ~("»,n)  o  0  (m,n)  ~(m,n) 


1  1 


T  ; 


61  =  yy  (-  — )6Y*"  ,  /  /  a  a*  6  f^-Od^dn  Y  . 

'  11  LL  T  ~(ni,n)  q  ;0  -  ~,n  (m,n)  ~(m,n) 


1  1 


61 


L  T 

T  -  1 


19  =  I  (f){<SY*  [k,  /  a(0,n)a  (0,n)dn 

12  n=l  L  ~(1’n)  0  ~ 


+  k2K2  J  a  ^ a  (0 ,ri)dr)] Y ^ j  ^ 


+  6Y*3  [k,  /  a(i,n)al(i»n)dn 

~(K,n)  3  o  ~ 


Also 


*  V2  /„  ;,5(1'n)Ac(1'n)dn^(K,n)1 


K  k 


61  =  l  —■  6Y *  ,  .  /  aK,l)aT(5,0)dC  Y 

13  m=l  K  ^m'LJ  0 

«!  -  I  I  (-  Ig6fa)5I^  nl  /  /  a(5,n)d5dn 

1  m=l  n=l  KL  (m>n)  0  0 


(5-11, part  B) 


|  E  C-  Tgmc— )6Y^  I1  J1  a(5,n )6fm  (5-5)d5dn 

2  m=l  n=l  L  0  0  ~  (m*n^ 


6J3  =  X  &>6Yf  t  i  ^  ^  Y  «)a(5fl)d5 
3  m=l  K  (m’L'  0  0  Cm).  ~ 


(5-12) 
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I 

With  61 £ ,  6j-  of  the  variational  statement  Eq.  (4-1)  written  in  terms  of 
a(£,q),  Y  ana  6Y*  as  given  in  Eqs,  (5-11)  and  (5-1G),  we  can  now  discuss 
the  construction  of  various  element  matrices  and  force  vectors. 

}  * 

5.2  Common  Element  Matrices.  A  common  element  matrix  is  one  which 
is  the  same  for  all  elements.  There  are  four  such  matrices  in  this  solu¬ 
tion  formulation  and  they  can  be  identified  in  the  expressions  for  61,, 
fil2,  <$I3  in  Eqs.  (5-11)  .  Let: 


K  L  »pv  s  t 

6Ii  =  I  l  T~  6Y,  ,  B  Y,  , 

1  m=l  n=l  L  ~Cm,n)  ~  ~(m,n) 


61,  = 


K  L 

l  l  (P-gsinot)  f  6Y  f  [(K-m+l)A-D]Y  ^ 

m=l  n=l  L  ~  ~ 


61-  »  £  l  (-  _L.)6Y*T  CY 

3  m=l  n=l  TK  ~Cra»n)  ~  ~Cm,n) 


(5-13) 


Then 


1  1 


A  = 


B  = 


/  /  3  r®  p  d^dq 

0  0 


1  1 


/  /  a  rra  rr  d£dq 

0  0 


1  1 


C  =  /  /  a  a1  d£dn 


0  0 


D  = 


11  t 

/  /  d£dq 


0  0 


(5-14) 


Each  of  these  matrices  A,  B,  C  or  D  is  of  size  16  x  16  and  thus  has  256 
values  to  be  evaluated.  Since  they  are  independent  of  the  element  para¬ 
meters,  they  need  only  to  be  computed  once  for  all  the  elements. 
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5.3  Element  Matrices  Due  to  Discontinuous  Load  Function.  These  ele¬ 
ment  matrices  are  identified  with  the  terms  of  6 I 4  through_6ljj  in  Eqs. 

(  r  ii\  .  .1 _  J  ‘  -  _  -  a.  .* I  f  ( r*  r  'i  J  /’C  _  r\  -**  fKni  r 


discontinuity  is  defined  by  £  =  £(r|).  The  specific  form  of  £(n)  is  now 
given.  Since  we  assume  that  the  projectile  moves  with  a  constant  accel 
eration  B  and  that  the  velocity  and  position  are  zero  at  t  =  0,  the  pro 
jectile  position  can  be  written  as 


x  =  x  =  —  Bt2 
2 


(5-15) 


From  Eqs.  (5-4),  one  has 


(5-16) 


Eq.  (5-16)  is  obviously  dependent  on  the  element  parameter  m,n.  In  addi¬ 
tion,  the  integrations  of  6I4  through  6 1 ^ ^ ,  will  depend  on  whether  and  how 
the  curve  of  Eq.  (5-15)  actually  goes  through  the  particular  element  (m,n). 
Thus  these  element  matrices  must  be  evaluated  element  by  element. 

5.4  Boundary  Element  Matrices  and  Non-Self- Ad joint  Matrices.  The 
contribution  due  to  the  end  elements  are  identified  with  the  <S  1^2  an^  ^*13 

terms.  Let 


Then,  it  is  clear  from  Eqs.  (5-12)  that 


5i  =  /  a(o,n)aT(o,n)dn 


0  ~ 


-1  T 

B  2  =  ]  a^(l,q)a,^(l,n)dr) 
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and 


1 

B,  =  /  a(l,n)a  (l,n)dn 
0  ~ 


T 

54  =  /  a^ci,n)a>c(i,n)dn 


B-  =  /  a(£,l)a(£,0)d£ 

0  ~ 


It  should  be  mentioned  that  I13  is  the  only  term  in  the  present  problem 
which  contributed  to  the  non-self-adjointness  of  the  problem.  This  is  also 
manifested  in  the  element  matrix  By  which  is  the  only  nonsymmetric  matrix. 

5.5  Element  Force  Vectors.  These  force  vectors  are  defined  in  con¬ 
junction  with  6Jj,  6J2  and  6J3  in  Eqs.  (5-12).  Let 

6J,  -  I  I  (-  Ig.c-P?0L)6Y*T  .  F, 

1  m,l  n=l  KL  ~<ra'n)  U 

K  L 


6J„ 


=  l  l  (-  IEI!1££§«)6y*T  f 


m=l  n=l 


~(m,n)  ~2 


K  k7  T 

*  J,  (T>  5* 


m= 


1  1 

h  ~  I  /  a(£,n)d£dn 


0  0 


1  1 

F7  =  /  /  a(S,n)<5 

0  0 


(m,n) 


(5-5)d?dn 


F3  ’  /  YW(5)a(5,l)dE 


5.6  Flow  Chart  of  the  Computer  Program.  As  shown  in  Figure  5  the 
computations  in  the  present  program  consist  of  two  parts:  data  prepara¬ 
tion  and  solution  program. 
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FIGURE  5.  Flow  Chart  of  the  Computer  Program. 


6.  DEMONSTRATIVE  RESULTS  OF  COMPUTATIONS.  Some  preliminary  results 
have  been  obtained  by  the  computer  program  described  in  previous  sections. 
Our  present  objective  is  simply  to  demonstrate  that  the  concept  of  the 
unconstrained  variational  formulation  works  for  the  initial  boundary  value 
problem  of  gun  dynamics,  that  the  initial  data  can  indeed  be  recovered  by 
the  unconstrained  process  and  that  the  nonconservative  and  highly  discon¬ 
tinuous  nature  of  the  loads  can  be  tested  in  a  routine  and  rather  elegant 
manner.  Thus  the  numerical  values  obtained  are  not  considered  final. 

Their  verifications  will  be  carefully  performed  and  presented  in  a  report 
to  follow. 

The  numerical  values  for  those  nondimensional  parameters  listed  in 
Section  2  are  given  here. 

T  =  0.07695 
P  =  4.8025 
m  =  0.012 
3  =  308.79 
a  =  0.5236 

g  =  0.020 

These  data  are  calculated  from  an  approximate  model  for  a  105  mm  M68  cannon 
with  the  following  material  and  dimensional  parameters: 

E  (Young's  modulus)  =  30xl06  psi 

p  (density)  =  0.2818  lb/in3 

I.D.  (of  gun  tube)  =  105  mm  =  4.13" 

O.D.  =  7.28"  (average  of  5.69"  muzzle 
O.D.  and  8.88",  breech  end  O.D.) 

A  =  28.23  in2 

I  =  123.60  in4 

l  =  17.5'  =  210" 

„  1/2 

C  =  r£~~_)  =  0.10396  sec. 

El 
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nip  (projectile  mass)  =  210  lb 

P  (bore  pressure)  =  30,000  psi 

Bp  (projectile  acceleration)  =  500,000  ft/sec2 

h  initial  position  of  the  tube  is  assumed  to  be  the  static  deflection  of  the 
tube  under  gravitational  force  q  =  pAg.  Hence 

Y (x)  =  -1-  (x4-4x3+6x2) 


where  q  =  and  g  is  the  gravitational  acceleration.  Hence  q  =  0.01987. 

Using  the  data  given  above,  the  solution  of  a  test  problem  is  shown  in 
Tables  hi  and jy.  The  exact  values  of  the  initial  data,  which  corresponds  to 
the  static  deflection  of  the  tube,  are  given  in  the  parentheses.  From  these 
tables,  it  is  clearly  shown  that  the  given  initial  data  have  indeed  been 
recovered  in  the  solutions. 

TABLE  hi.  CALCULATED  y (x, t) (xlO' 3)  FROM  THE  PRESENT  COMPUTER  PROGRAM 
(EXACT  INITIAL  DATA  IN  PARENTHESES) 


f  K 

0 

0  2 

0.4 

0.6 

0.8 

0.10 

0 

0.00000 

0.17477 

0.60803 

1.18801 

1.83469 

2.49996 

(0.00000) 

(0.17467) 

(0.60800) 

(1.18800) 

(1.83467) 

(2.50000) 

0. 5T 

0.00000 

0.1S504 

0.59019 

1.17861 

1.82291 

2.48327 

i 

0.00000 

0.21705 

0.58048 

1.15876 

1.77023 

2.45612 

TABLE  IV.  CALCULATED  (xlO-3)  FROM  THE  PRESENT  COMPUTER  PROGRAM 

dx 

(EXACT  INITIAL  DATA  IN  PARENTHESES) 


0 

0.2 

0.4 

0.6 

0.8 

0.10 

0 

0.00000 

1.62731 

2.61662 

3.12140 

3.30629 

3.32996 

(0.00000) 

(1.62667) 

(2.61333) 

(3.12000) 

(3.30667) 

(3.33333) 

0.5T 

0.00000 

1.63528 

2.62714 

3. 14336 

3.27633 

3.31506 

T 

0.00000 

2.06887 

1.68236 

3.03449 

3.23987 

3.55142 
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COMMANDER 

US  ARMY  HARRY  DIAMOND  LAB 
ATTN:  TECH  LIB  1 

2800  POWDER  MILL  ROAD 
ADELPHIA,  MD  20783 

DIRECTOR 

US  ARMY  INDUSTRIAL  BASE  ENG  ACT 
ATTN:  DRXFE-MT  1 

ROCK  ISLAND,  H  61201 

CHIEF,  MATERIALS  BRANCH 

US  ARMY  R&S  GROUP,  EUR  1 

BOX  65,  FPO  N.Y.  09510 

COMMANDER 

NAVAL  SURFACE  WEAPONS  CEN 

ATM:  CHIEF,  MAT  SCIENCE  DIV  1 

DAHLGREN,  VA  221,48 

DIRECTOR 

US  NAVAL  RESEARCH  LAB 
ATM:  DIR,  MECH  DIV  1 

CODE  26-27  (DOC  LIB)  1 

WASHINGTON,  D.C.  20375 

NASA  SCIENTIFIC  &  TECH  INFO  FAC 
P.O.  BOX  8757,  ATM:  ACQ  BR  1 

BALTIMORE/WASHINGTON  INTL  AIRPORT 
MARYLAND  21240 


NO.  OF 
COPIES 

COMMANDER 
DEFENSE  DOCU  CEN 

ATTN:  DDC-TCA  12 

CAMERON  STATION 
ALEXANDRIA,  VA  22314 

METALS  &  CERAMICS  INFO  CEN 
BATTELLE  COLUMBUS  LAB  1 

505  KING  AVE 
COLUMBUS,  OHIO  43201 

MPDC  1 

13919  W.  BAY  SHORE  DR. 

TRAVERSE  CITY,  MI  49684 

MATERIEL  SYSTEMS  ANALYSIS  AC TV 
ATM:  DRXSY-MP  1 

ABERDEEN  PROVING  GROUND 
MARYLAND  21005 


NOTE:  PLEASE  NOTIFY  COMMANDER,  ARRADCOM,  ATTN:  BENET  WEAPONS  LABORATORY, 
DRDAR— LCB— TL ,  WATERVUET  ARSENAL,  WATERVLIET,  N.Y.  12189,  OF  ANY 
REQUIRED  CHANGES. 


EXTERNAL  DISTRIBUTION  LIST 


NO.  OF 
COPIES 

ASST  SEC  OF  THE  ARM! 

RESEARCH  &  DEVELOPMENT 
ATTN:  DEP  FOR  SCI  &  TECH 
THE  PENTAGON 
WASHINGTON,  D.C.  20315 

COMMANDER 

US  ARM!  MAT  DEV  &  READ.  COMD 
ATTN:  DRCDE  1 

5001  EISENHOWER  AVE 
ALEXANDRIA,  VA  22333 

COMMANDER 
US  ARM!  ARRADCOM 

ATTN:  DRDAR-TSS  2 

DRDAR-LCA  (PLASTICS  TECH  1 
EVAL  CEN) 

DOVER,  NJ  07801 

COMMANDER 
US  ARMY  ARRCOM 

ATTN:  DRSAR-LEP-L  1 

ROCK  ISLAND  ARSENAL 
ROCK  ISLAND,  LL  61299 

DIRECTOR 

US  ARMY  BALLISTIC  RESEARCH  LABORATORY 
ATTN:  DRDAR-TSB-S  (STINFO)  1 

ABERDEEN  PROVING  GROUND,  MD  21005 

1 

COMMANDER 

US  ARMY  ELECTRONICS  COMD 
ATTN:  TECH  LIB  1 

FT  MONMOUTH,  NJ  07703 

COMMANDER 

US  ARMY  MOBILITY  EQUIP  R&D  COMD 
ATTN:  TECH  LIB  1 

FT  BELVOIR,  VA  22060 


NO.  OF 
COPIES 

COMMANDER 


WARREN,  MICHIGAN  48090 
COMMANDER 

US  MILITARY  ACADEMY 

ATTN:  CHMN,  MECH  ENGR  DEPT  1 

WEST  POINT,  NY  10996 

COMMANDER 
REDSTONE  ARSENAL 

ATTN:  DRSMI-RB  ? 

DRSMI-RRS  1 

DRSMI-RSM  1 

ALABAMA  35809 

COMMANDER 

ROCK  ISLAND  ARSENAL 

ATTN:  SARRI-ENM  (MAT  SCI  DIV)  1 

ROCK  ISLAND,  IL  61202 

COMMANDER 

HQ,  US  ARMY  AVN  SCH 

ATTN:  OFC  OF  THE  LIBRARIAN  1 

FT  RUCKER,  ALABAMA  36362 


COMMANDER 

US  ARMY  FQtf  SCIENCE  &  TECH  CEN 
ATTN:  DRXST-SD  1 

220  7TH  STREET,  N.E. 

CHARLOTTESVILLE,  VA  22901 

COMMANDER 

US  ARMY  MATERIALS  &  MECHANICS 
RESEARCH  CENTER 

ATTN:  TECH  LIB  -DRXMR-PL  2 

WATERTOWN,  MASS  02172 


US  ARMY  TANK-AUTMV  R&D  COMD 
1  ATTN:  TECH  LIB  -  DRDTA-UL  1 

MAT  LAB  -  DRDTA— RK  1 


NOTE:  PLEASE  NOTIFY  COMMANDER,  ARRADCOM,  ATTN:  BENET  WEAPONS  LABORATORY, 
DRDAR-IGB-TL,  WATERVLIET  ARSENAL,  WATERVLIET,  N.Y.  12189,  OF  ANY 
REQUIRED  CHANCES. 


